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Abstract 

This article describes a complex network model 
whose weights are proportional to the difference 
between uniformly distributed "fitness" values as- 
signed to the nodes. It is shown both analytically 
and experimentally that the strength density (i.e. 
the weighted node degree) for this model, called 
derivative complex networks, follows a power law 
with exponent 7 < 1 if the fitness has an upper 
limit and 7 > 1 if the fitness has no upper limit 
but a positive lower limit. Possible implications for 
neuronal networks topology and dynamics are also 
discussed. 

1 Introduction 

Great part of the interest focused on complex net- 
works ^ |2 El recently stems from scale free or 
power law distributions of respective topological 
measurements, such as the node degree. At the 
same time, weighted complex networks have at- 
tracted growing interest because of their relevance 
as models of several natural phenomena, with spe- 
cial attention given to systems used for distribu- 
tion/collection of materials or information. 

There are two main ways to approach the de- 
grees of weighted networks: (i) by thresholding the 
weights and using the traditional node degree pQ; 
and (ii) by adding the weights of the edges attached 
to each node, yielding the respective node strength 
Related recent works include the identification 
of strength power law in word association networks 
[S] studies of scientific collaborations and air- 



transportation networks |7j, the analysis of amino 
acid sequences in terms of weighted networks [§], 
the investigation of weighted networks defined by 
dynamical coupling between topology and weights 
0] , analytical characterization of thresholded net- 
works [B], as well as the characterization of motifs 
|Hj and shortest paths in weighted networks 

The current article describes a network whose 
nodes have a respective "fitness" value and every 
node is connected to all other nodes with higher 
fitness through an edge whose weight is determined 
by the difference between the fitness of the nodes. 
Because of such a construction, these networks are 
henceforth called derivative complex networks. Pre- 
vious works considering the difference of fitness val- 
ues associated to the network nodes include the 
image analysis approach reported in |llj , which in- 
volves the thresholding of the fitness difference val- 
ues in order to perform image segmentation and 
the gradient networks discussed in |1'2I H3) , which 
take into account the edges corresponding to the 
highest differences, i.e. to the gradient of a scalar 
field associated along the network nodes. 

The developments in the current article represent 
a continuation and extension of preliminary investi- 
gations reported in We show both analytically 
and numerically that the strength densities for the 
derivative network model follow a power law with 
exponent 7 < 1 if the fitness has an upper limit 
and 7 > 1 if the fitness has no upper limit but a 
positive lower limit. 

We start by describing the derivative network 
model and follow by calculating its respective 
strength distribution and discussing possible impli- 
cations for neuronal networks. 



2 Derivative Networks loss of generality, to be unitary: < ip < 1. For a 

power law distribution 

Consider a network with N nodes. To each node i, 

a "fitness" value (pi > is randomly assigned fol- p(ip\ = cnp a ~ x (6) 

lowing a distribution p{ip). The fitness derivative 

determines the connectivity of the nodes in the net- to be a valid probability distribution in the interval 
work. A directed arc linking node j to node i is < ip < 1, we must have a > 0. From Eqs. J2J|, © 
drawn iff ipi > tpj; in that case, the arc has weight and © we have: 
given by a function of the fitness difference: 

afa - (1) S = N ^ aar^Ktp xf^ dx, (7) 
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We are here interested in the study of the distribu- resulting in: 
tion of input strengths in the network under some 

assumptions on the functions p and a. The input s = NbaB(a, f3)(p a+ ^^ 1 , (8) 

strength of node i, Si, is the sum of the weights of 

the arcs linked to i: where B(a, b) is the Euler beta function: 

Si=^2wji, (2) 



B(o,6) = / t a - 1 (l-t) b - 1 dt. 
Jo 

Note that this expression implies an upper limit for 

From Eqs. (0) and QJ, the strength of a node is S ^ ven by 
determined by its fitness ip: Smax = NbaB(a,(3). (9) 

, \ i is /o\ The probability density of s, p(s), is given by: 

= a i<p - <p ); ( 3 ) 



where the sum is taken over all nodes j that have 
an arc linked to node i 



¥'<<p 



P(*) = P(<P)% (10) 
where the sum is taken for all nodes with ip' smaller 

than ip. Considering the distribution of fitness r^i j ■ , ■ « -uj. • j r n m . • 

j; . b 1 he derivative can be obtained from Eq. giving: 

/5(</?), for the limit of large number of nodes N we 
can write: 



s = N I p(x)a(ip - x) dx. (4) 



This expression is only valid for < s < s max , and 
For the weight function, we assume a power law is a power law p(s) ~ s~ 7c/ with 



as: 



fffa:) = fox' 3 - 1 . (5) P- 1 Mo v 

v ; w 7c/ = — r-3 — 7- ( 12 

a + p — 1 

We are interested in the case where the weight 

grows with the difference in fitness, implying (3 > 1. Considering a > and f3 > 1 we have < 7 < 1. 

For the fitness distribution, a power law distribu- An interesting special case is that of uniform d is- 

tion will be also considered. In the following sec- tributianj which corresponds to a = 1 in Eq. ®, 

tions, distributions with an upper limit and with a resulting in 
positive lower limit will be considered. ^ 

2.1 Upper bound fitness 

In case we also have that the weight is proportional 

First we study the case where the values of fitness to the fitness difference, i.e. (3 — 2 in Eq. JSJ, then 

have a maximum value that we assume, without ju = h. 
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2.2 Lower bound fitness 

We consider now the case where the fitness value 
has no upper limit, but a positive lower limit, that 
we assume, without loss of generality, to be unitary. 
In this case, we write the probability density as 

Krf^^ 1 , (13) 

which is correctly normalized in the interval 1 < 
<p < oo for a > 0. 

Using Eqs. ®, 0, and {I3J we have: 

ax- a - x b(<p - xf- 1 dx, (14) 

giving: 

s = JV6aB(-,l,-a,/9)^~ 1-a . (15) 

<fi 

where B(x, y, a, b) is the generalized incomplete 
beta function: 

rv 

B(x, y, a, b) = / i 0_1 (l - 1) 6 " 1 dt. 

J X 

For ip w 1, B(i, 1, — a,/?) « giving s « 0. For 
large <p, B(~, 1, —a, /3) ~ (p a . Then 

iVi»aB(-,l,-Q!,j9) ~c</? a (16) 

and therefore we get: 

s ~ c^ 1 (17) 

for large 95. 

Proceeding in a similar fashion to the previous 
case we get: 

P( s ) - a 7 c ^ _ls 3=1 > ( 18 ) 

P - 1 

valid for large s. 

Considering that j3 > 1 and a > 0, this expres- 
sion corresponds to a power law p(s) ~ s -71 with 
Tl > !• 

3 Simulation results 

We simulated derivative networks with upper and 
lower limit for three sets of a and parameters 
each. Figure ^ shows the strength distribution for 



networks with upper fitness limit; figure [3 shows 
the results for networks with lower limited fit- 
ness. The results assumed N = 2000 nodes and 
b = 0.0005; a total of 100 networks were considered 
for each set of parameters, and the results show the 
mean values and respective standard deviations. 

The curves obtained by the simulations resulted 
remarkably close to the analytical results for upper 
limited networks, except for the smallest strengths. 
This is because the strength distribution p(s) ~ 
s 1 approaches a discontinuity at zero. 

Less adherence between theoretical and experi- 
mental values was observed for the lower limited 
networks as a consequence of the fact that the ap- 
proximation in Eq. (|16fl is only valid for large val- 
ues of s, which corresponds to small values of the 
probability p(s), resulting in poor statistics. 



4 Discussion 

We have described two simple models of derivative 
networks and shown that their node instrength dis- 
tributions follow a power law under some assump- 
tions on distribution of fitness and weight assign- 
ment. This result implies that although several 
nodes of a derivative network have low strength 
values, there are hubs characterized by a variety 
of large instrength values. Because the instrength 
of a node can be understood as the total weighted 
influence it receives from the adjacent nodes, it fol- 
lows that in case the node state is directly pro- 
portional to its instrength, the linear dynamics of 
such networks will be characterized by a near power 
law with similar parameters as those of the in- 
strength density. Such an interpretation is par- 
ticularly interesting from the perspective of inte- 
grating neuronal network and complex network re- 
search (e.g. 2HDi where each neuron is represented 
by a node and the synaptic connections by directed 
edges whose weights reflect the respective synaptic 
strengths. The "fitness" values in this case could 
be related to gradient of neurotrophic growth fac- 
tors or depolarization bias facilitating the action 
potential |161 IT?] . Another approach worth further 
attention would be to consider derivative networks 
as models of perceptual and cognitive processes, 
which are known to explore derivatives as a means 
of diminishing stimuli redundancy (e.g. |18j'). 
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Figure 1: Strength distribution in upper limited 
networks. Experimental results are represented by 
points with error bars (one standard deviation) and 
the continuous curves are the theoretical expres- 
sions. Results for: (a) a = 2, (3 = 2 (7 = 1/3); 
(b) a = 1,(3 = 2 (7 = 1/2); (c) a = 1,(3 = 3 
(7 = 2/3). 



Figure 2: Strength distribution in lower limited 
networks. Experimental results (points with error 
bars) are represented together with a continuous 
line displaying the theoretical asymptotic inclina- 
tion. Results for: (a) a = 0.5,(3 = 1.5 (7 = 2); 
(b) a = 1.5,13 = 2 (7 = 2.5); (c) a = 1,(3 = 1.5 
(7 = 3). 
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